Chiral anomaly and strength of the electron-electron interaction in graphene 
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The long standing controversy concerning the effect of electron - electron interaction on the 
electrical conductivity of an ideal graphene sheet is settled. Performing the calculation directly in 
the tight binding approach without the usual prior reduction to the massless Dirac (Weyl) theory, 
it is found that, to leading order in the interaction strength a — e 2 /hvo, the DC conductivity 
a /(To = 1 + Col + O (a 2 ) is significantly enhanced with respect to the independent-electrons result 
cro, i.e. with the value C = 0.26. The ambiguity characterizing the various existing approaches is 
nontrivial and related to the chiral anomaly in the system. In order to separate the energy scales in a 
model with massless fermions, contributions from regions of the Brillouin zone away from the Dirac 
points have to be accounted for. Experimental consequences of the relatively strong interaction 
effect are briefly discussed. 
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Introduction. It has been demonstrated recently that a 
graphene sheet, especially one suspended on leads, is one 
of the purest electronic systems. The scattering of charge 
carriers in suspended graphene samples of submicron 
length is so negligible that the transport is ballistic 0, ■ 
The novelty of the physics of undoped graphene is in the 
ability to probe the "ultrarelativistic" physics of excita- 
tions leading to numerous similarities with phenomena 
previously associated with the high energy physics. Ex- 
amples include Zitterbewegung and Klein tunneling 0], 
electron - hole (Schwingcr) pair creation by an electric 
field Q, a possibility of dynamical (chiral) symmetry 
breaking by electron interaction effects [j| (exciton con- 
densation) and the chiral (parity) anomaly Q. The lat- 
ter, a quantum anomaly, attributed to graphene long be- 
fore its discovery, is one of the most remarkable features 
of a relativistic field theory with massless fermions 0, Q. 
Generally it is associated with the fact that a classical 
symmetry is "broken" by quantization in the case of an 
infinite number of degrees of freedom, when the ultravio- 
let (UV) cutoff is necessary. Chiral anomaly means that 
the classical axial U (1) symmetry is violated. This led to 
explanations of such physical phenomena 0] as 7r° — > 77 
decay (that would be suppressed by the symmetry), the 
solution of the problem of the large mass of the 77 meson 
(zero, if it were to be a Goldstone boson of a nonanoma- 
lous symmetry) etc. The anomalies are notorious in that 
calculations of a well defined physical quantity using dif- 
ferent UV cutoff procedures (for example the sharp mo- 
mentum cutoff, lattice regularization or a properly de- 
fined dimensional regularization) led to different finite 
values. The physical essence of this ambiguity is that 
there is no simple separation between the UV and in- 
frared (IR) physics and certain care should be exercised 
in construction of the correct effective low energy model. 
This might be suspected to occur in theory of graphene. 
In description of graphene, while the starting point might 



be an atomic or tight binding modcl|9(, one typically 
replaces it by an massless effective Dirac (Weyl) model 
"near" its two Dirac points constituting the Fermi "sur- 
face" of undoped graphene. 

In this note we point out that the elucidation of the am- 
biguities encountered in the theory of the (apriori strong) 
Coulomb interactions should be associated with a careful 
treatment of the separation of scales due to the anomaly. 
We show in detail, using the tight binding model pro- 
viding a natural UV cutoff, that some aspects of the 
graphene physics are not dominated by the two Dirac 
points of the Brillouin zone at which the spectrum is 
gaplcss. The low frequency conductivity in the quasi- 
dielectricphase below the exciton condensation critical 
coupling [E EH a = e 2 /hv < a c (neglecting weak loga- 
rithmic renormalization of the electron velocity [1CH13| , 
v = vq ~ 10 6 m/s), is given in terms of its value in the 
noninteracting theory, ctq = e 2 /Ah, by 



a (uj) /a = l + Ca + O (a 2 ) 



(1) 



This expression is valid for frequencies below the hopping 
energy 7 = 2.7eV. The static dielectric constant is there- 
fore given by e" 1 = 1 - na/2 (1 + Ca) +O (a 3 ). The 
value of the only numerical constant C appearing here 
has been a matter of intense controversy. The first de- 
tailed calculation by Herbut, Juricic and Vafckjl4| utiliz- 
ing a sharp momentum cutoff regularization of the Dirac 
model provided a value of order 1: 
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0.51. 
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The use of the sharp momentum cutoff was criticized by 
Mishchenko Il5j , who obtained a value of 
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making a "soft" momentum cutoff regularization. He 
supported this choice by the consistency of the Kubo 
and the kinetic equation calculations of conductivity with 
that of the polarization function (dielectric constant). 
The consistency required a modification of the long range 
interaction so that it becomes UV cutoff dependent. It 
was further supported by Sheehy and Schmalianflij who 
used yet a different cutoff procedure and pointed out that 
only such a small value of C can explain the experimen- 
tal observation of the optical conductivity in graphene 
on a substrate 17] , which is within 1% of ctq- This ap- 
parently closed the issue. Albeit such a small numerical 
value would have profound physical consequences even 
beyond the transport and dielectric properties. 

Nevertheless the interaction strength C was recalcu- 
lated once again by Vafek, Juricic and Herbut[l8| who 
argued that the modification of the interaction requires 
simultaneously a Pauli - Villars regularization of mass- 
less fermions. They applied yet another regularization, 
making the space dimensionality fractional, D = 2 — e 
(similar to the space - time 4 — e regularizations that has 
been long in use in high energy and critical phenomena 
physics [8() that modified both the current operator and 
the interaction in such a way that they satisfy the Ward 
identities and obtained 
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0.26. 
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The dimensional regularization is questionable on physi- 
cal grounds and in a comprehensive subsequent work[19j 
the authors reaffirmed the small value and it seems 
that it is a commonly accepted one. To refute the ear- 
lier calculation of ref.[l8| they write "of course, satisfy- 
ing the Ward-Takahashi identity does not guarantee that 
the regularization scheme will produce the exact value of 
C for the physical system. We believe that if a really 
quantitative result is desired for the constant C, then 
one should resort to a complete electronic structure cal- 
culation (based, for example, on a realistic tight-binding 
Hamiltonian) rather than working with an effective low- 
energy theory". We followed this path, but surprisingly 
found that the tight binding value is . The situation 
is further complicated by other values in literature like 
C = 0.34 obtained in a dielectric constant calculation (20|. 

To reveal the origin of the ambiguity exhibited by the 
various values of C (there is a consensus that all the 
calculations are mathematically sound [3]), we use a dy- 
namical approach (used previously to address the ambi- 
guity of the noninteracting case 

dEH) directly in the DC 
case by " switching on" a uniform electric field in the tight 
binding model with Coulomb interactions, and then con- 
sidering the large-time limit. This approach (known in 
field theory as the " infinite hotel story" ) is the best way 
to reveal physical effects of anomalies 0, Si- O ne can di- 
rectly separate the contributions from the neighborhood 
of Dirac points and the " anomalous" contributions from 
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FIG. 1: Honeycomb lattice for graphene. The sublattice A 
(red) is spanned by the lattice vectors ai,2 = § (±1, V3) 
where a ~ 3A. The three nearest neighbours on sublat- 
tice B (blue) are displaced by di = i (ai — ag), = 
| (ai + 2a2), d3 = — | (2ai +a2). Wilson links W describ- 
ing the minimal coupling to the vector potential A (r, t) are 
defined in Eq.©. 



the rest of the Brillouin zone, so that one can decide what 
regularization of the effective Weyl theory is the correct 
one. We have also performed a standard diagrammatic 
Kubo formula calculation of the general AC conductivity 
within the tight binding model and obtained the same 
result. 

The tight binding model and its linear response to an 
electric field. Electrons in graphene are described suffi- 
ciently accurately for our purposes by the 2D tight bind- 
ing model of nearest neighbour interactions in external 
field described by Wilson links 22j : 



K[A] = -~/^LdgWn,i<Z B+dt + hci 



(5) 
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Here A (r, t) is the vector potential and dp creates an 
electron with spin a (summation over a implied) on the 
sites of the honeycomb lattice r n = riiai + n2&2, where 
lattice vectors ai^ and the nearest neighbours displace- 
ments di are defined in Fig.l. 

Coulomb interactions between electrons are 
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'jr. The corresponding current density operator (in 
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Heisenberg picture) is cSK [A] /SA (r, t): 

J(r,t) = -^Vd t /" Sfr-i^-sd;) (7) 
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This describes a network-like flow of currents on links be- 
tween neighboring sites in Fig.l. As was emphasized in 
the context of quasi-local interaction models (Hubbard 
models) in graphene in ref. [231] (and much earlier in the 
context of lattice gauge models of particle physics (22^) 
this model satisfies all the (nonanomalous) Ward iden- 
tities associated with charge conservation and therefore 
no trial-and-error modification of the current operator is 
needed. 

Let us consider a uniform electric field along the y 
direction E = — J- (0, A (t) /c) switched on at t = 0. 
The current density is expanded to first order in A as 
J = + J p , with the relevant components being, 
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Averaging the expectation value of current density over 
the sample area S, j (t) = ^ J r (<p \J V (t) \ </>), one obtains: 
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The time independent Hamiltonian H = K + V, K = 
K [A = 0] and its ground state \<f>) are expanded to first 
order in the interaction V. The tight binding model K 
has a spectrum = ± |/ik| determined by the structure 
function of the links hy = -7^, e~ lk d \ The DC field is 
defined by A (t) = —cEt and resuls of direct calculation 
are presented and discussed in what follows. 

The evolution of the current. The current density to 
first order in interactions is 



j (t) = o a E [C° (t) + aC (t)] 



(10) 



The components of the minimal dimensionless conduc- 
tivity C° (t) are written as integrals over the Brillouin 
zone (BZ), see Fig. 2: 




FIG. 2: Contours of the tight binding energy £k in the entire 
Brillouin zone of the honeycomb lattice. Circles of radius K, 
around the two Dirac points are the parts described by the 
effctive low energy Weyl model. 
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where (k = In (z^h'^) , z^ = h^/s^, and a prime denotes 
a derivative with respect to momentum along the field, 
k y . 

The conductivity includes two apparently linearly di- 
vergent in time "acceleration" parts. However their sum 
integrated over the BZ vanishes since it is a full deriva- 
tive of a periodic function, Re ^ k (h^h^ + /e^ = 
4^] k ej| = 0. This cancellation, albeit, is nontrivial: con- 
tributions come from the whole BZ. When one uses the 
effective Dirac theory, h^ = voh(k x ± ik y ) = vohke ±lLPk , 
Ck ~ foftcos</?k and integration over each of the two cir- 
cles (with the cutoff radius /C) in Fig. 2 gives a positive 
UV divergent result: 
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ip = nvohK,. 
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This is canceled exactly by contributions from regions of 
BZ far from Dirac points in which the low energy effective 
model is not valid[2J|. Now that the "acceleration" parts 
have cancelled, the oscillating term in j p at large t limit 
can be safely calculated from the effective low energy 
theory. Indeed averaging over long times T = l/r/ by 
C = rj J t °° C (t) e~ vt , one obtains 
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as expected. Now we turn to the interaction corrections. 

Similarly as before, the linear in t "acceleration" cor- 
rections, 
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cancel each other beyond the Weyl model applicability 
domain, leaving oscillating terms of C p that now take a 
form (averaged over large times): 
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The integrals arc again computed (see Supplemental ma- 
terial for details) using the Dirac point approximation. 
This is the main result of the present work. In this man- 
ner we also calculated the AC conductivity of the tight 
binding model and results will be presented elsewhere. 

Summary and discussion. To summarize, we have cal- 
culated the electron-electron interaction contribution to 
DC and AC conductivity of undoped graphene within 
the tight binding model. Thus the controversy of what 
is the actual magnitude (even order of magnitude) of the 
corrections is resolved in favour of the intermediate value 
of the constant C = C^ 3 ) >> C (2 ). It is shown that the 
ambiguity between the three values originates in a non- 
trivial feature of massless fermions, the chiral anomaly. 
The major complication that massless fermions cause is 
the absence of a perfect scale separation between high 
energies (on atomic scale 7) and low energies (effective 
Weyl theory on the condensed matter scale << 7). We 
demonstrated that some aspects of the linear response 
physics are not dominated by the two Dirac points of the 
Brillouin zone at which the effective low energy model 
is valid. For example, large contributions (infinite, when 
the size of the Brillouin zone is being considered infinite) 
to the conductivity from the vicinity of the Dirac points 
arc cancelled by contributions from the region between 
them. Another famous consequence of this scale nonsep- 
aration is the "species doubling" of lattice fermions [22[, 
which in the context of graphene means that there nec- 
essarily appears a pair of Dirac points of opposite chiral- 
ity. The UV regularization of the effective theory does 
matter and, if one were to use such a model, the only 
regularization known to date to be consistent with the 
tight binding is the space dimensional regularization de- 
veloped in ref. [18| . The reason is not clear to us (espe- 
cially due to the fact that fully relativistic dimensional 
regularization in 2+1 anomalous theories is known to be 



problematic[i|,[22j]), but experiece with field theory would 
indicate that one can also construct a successful suffi- 
ciently simple Pauli - Villars kind of regularization. 

If the result C = were the correct one, the physics 
would look very different. Indeed such a small value 
would easily explain the experimenal absence of interac- 
tion correction s 1 16l in the AC conductivity of graphene 
on a substrate [l7|. The explanation probably resides 
elsewhere, for example in the dielectric constant of the 
substrate, screening due to puddles, etcjl^]. We have 
calculated the effect of screened interactions represented 
by the Hubbard model with quasilocal interactions (up 
to several nearest neighbours) and obtained a vanishing 
first order correction to AC conductivity at all frequen- 
cies in accordance with a general theorem [23j]. For local 
interactions this has been already noted in ref. [l9j ]. 

The intermediate value of C can also have a bearing on 
the putative exciton condensation due to strong Coulomb 
interaction that has not been yet experimentally observed 
even in suspended graphene samples Q and on interac- 
tion corrections to the dispersion relation of the exci- 
tations. The random phase approximation (RPA) and 
various large Nf results should be also derived from 
the tight binding model or from a properly regularized 
effective low energy one. It is well known in field the- 
ory that generally chiral anomaly effects appear only in 
one loop calculations Q and higher orders resummed 
in RPA or 1/Nf approximation should not lead to fur- 
ther ambiguities. It is remarkable to note that differ- 
ences between the values of C in equations Eqs. (|2|3|4j) 
are C*- 1,2 ) = ± 1/4. Sometimes due to anomalies 
similar differences in regularizations are related to cer- 
tain "topological invariant" stemming from the measure 
of the path integration over fermionic fields[§]. Here the 
situation is more complicated since we are dealing with 
correction due to interactions, not with the simple bubble 
diagram. 
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